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O^ ' The question of whether the pseudogap in high Tc cuprates is related to superconducting precursor 

effects or to the existence of extrinsic bosonic massive excitations is investigated on the basis of the 
Boson- Fermion model. The characteristic three peak structure of the electronic spectral function and 
the temperature dependent Fermi vector derived here are signatures for a two component scenario 
which can be tested by ARPES and BIS experiments. 



00 



o 



> 

O 
(N 

O 

a^ 



o 



X 
S3 



PACS numbers; 79.60.-i, 74.25.-q, 74.72.-h 



I. INTRODUCTION 



The recent direct experimental verifications Ep] of a pseudogap in the density of states of the electrons in the 
metallic layers of the cuprate high temperature superconductors in the underdoped regime is presently considered 
to be one of the key elements to guide us in unravelling the underlying mechanism for superconductivity in these 
materials. 

The most immediate temptation to explain the pseudogap which opens at a certain temperature T* , sometimes 
well above the superconducting transition temperature, is to consider it as a consequence of superconducting phase 
fluctuations corresponding to a finite and relatively long lived amplitude of the order parameter above Tc- Support in 
favor of such a scenario comes from the fact that the symmetry of the pseudogap is the same as that of the true gap 
and that the materials are quasi 2D structures, for which Cooper pairs are expected to exist in form of propagating 
states above Tc Isl. It is in this context that many theoretical papers base their study on the intricate interplay 
between electronic quasi-particles and their bosonic resonant states. 

Scenarios for the pseudogap based on systems with attractive interaction between the electrons exhibit such a pre- 
cursor pseudogap in the cross-over regime between BCS an a superfiuid state of preformed tighly bound electron pairs. 
Quite generally and independent of any particular scenario T^ is then determined by a temperature Tq characterizing 
the onset of phase fluctuations of the order parameter, while the meanfield critical temperature Tmf describes the 
formation of pairs. This leads to a suppression of intensity of low lying single particle excitations and thus to the 
appearance of a pseudogap Q . Specific scenarios attempting to describe such behavior are based on generic models 
for superconductivity with some effective attractive interaction [g-0 and on the negative U Hubbard problem gj^. 

Contrary to such scenarios implicating a direct precursor effect, others have been proposed where the pseudogap 
is induced by antiferromagnetic spin fluctuations rather than pair fluctuations. For single component systems such 
scenarios are based on strong repulsive interactions between the electrons. The pseudogap is then related to a SDW 
and/or Hubbard pseudogap close to a corresponding transition [ [lO||llf| . More precicely, the pseudogap features in 
the charge sector are derived from the pseudogap features of bosonic resonant states in the spin sector |l2-n3. A 
somewhat different scenario for the pseudogap is presented for the RVB physics, according to which chargeless spinous 
pair up at T* and holons condense at T^ p5[ . Similar in spirit are the SU{2) symmetry meanfield studies of slave 
boson approaches in strongly correlated systems [ [L6[ , which relate the pseudogap to the transition between small and 
large Fermi surfaces. Van Hove singularities also may play a role and were shown to stabilize d wave RVB states with 
pseudogaps near those singularities jlj. 

Apart from those single component scenarios, scenarios involving two component systems have been proposed. 
Experiments on the optical conductivity have for a long time suggested that in high Tc cuprates two types of carriers 



may be at stake, almost localized ones and itinerant ones |18|. Recent experiments on the specific heat are in support 



of such a hypothesis |19| and chemical analysis |20| favors the picture of metallic layers imbedded in highly polarizable 
dielectrics which could be the seat of localized polarons or bi-polarons. These results clearly point in the direction of 
two component scenarios such as described by the generic Boson-Fermion (B-F) model introduced by us along time 
ago |gl[ . On the basis of this model we have initially shown the opening of a pseudogap in the density of states of the 
Fermionic subsystem setting in below a certain temperature T* , the breakdown of Fermi liquid properties [E2U23] and 
the manifestations of that in the thermodynamic, magnetic and transport properties p4,M. A scenario related to 
that has been proposed on the basis of stripes in the metallic layers, where ID spin stripes alternate with ID metallic 
stripes [Eq| . The spin stripes produce bosonic entities in form of tightly bound spin singlets which are exchanged by 
pair hopping processes with electron pairs in the neighboring metallic stripes. 

The purpose of this communication is to elucidate the nature of the pseudogap as given in the scenario of the 
Boson-Fermion model. Is it due essentially to a precursor effect in low dimensional systems PTj or is related to a two 
particle resonant state triggered by the existence of extrinsic (in contrast to intrinsic in single component systems) 
bosonic modes such as localized bi-polarons which form part of the charge carriers of the system? For this purpose 
we shall try to interpret our previous numerical calculations [E2U23| in terms of an approximative analytic approach. 
This will permit us to intuitively understand the three peak structure for the spectral properties of the electronic 
subsystem, the reasons for the breakdown of Fermi liquid properties and indicate to what extent the opening of the 
pseudogap is related to a precursor to superconductivity. Moreover we shall show that for this model the Fermi vector 
will change as a function of temperature. This should in principle be measureable and thus help to clarify whether 
high Tc superconductors are one or two component systems. 



II. AN ANALYTIC APPROACH TO UNDERSTAND THE PSEUDOGAP 

The B-F model describes a system of localized Bosons hybridized with conduction electrons. The Bosons can be 
thought of as bipolarons (localized electron pairs with opposite spin) mainly located in the dielectric layers surrounding 
the Cu02 planes which contain the Fermions (conduction electrons) [ p8| . The corresponding Hamiltonian for this 
model is given by 

if =(zt-^)^C+C,;,-t J2 dc,a + (Ai3 - 2/.) 

X Y^ b+bi + vY^[btc^lC^^ + c+c+h] (1) 

i i 

where c\^ refers to the Fermion operators of the itinerant electrons and b] to the Boson operators. The spin 
indices are given by a, the bare hopping integral for the electrons is given by t, the Boson energy level by Ab and 
the Boson-Fermion pair exchange coupling constant by v. i denotes some effective site involving adjacent molecular 
clusters of the metallic and dielectric planes, which, as far as the metallic planes are concerned, should involve a unit 
of two neighboring CuO^ clusters and as far as the dielectric planes are concerned could be thought of two adjacent 
apex oxygens together with their nearby dopant atoms. Such pairs of electrons, corresponding to the term cXcl, in 
the metallic layers should be considered as intersite rather than onsite pairs. Correlation effects which are certainly 
important are totally neglected here, because we want to focus on physics introcuccd by pair hopping between the 
metallic and the dielectric layers. Since we are dealing here with a two-component system we have to impose a 
common chemical potential fi in order to guarantee charge conservation. 

We have previously evaluated the Fermion spectral functions and the density of states for this model for a one 
and respectively two dimensional bare (i.e., in the absence of Boson Fermion hybridization, v = 0) density of states 
using a selfconsistent lowest order fully conserving diagramatic procedure |g^,^. In the present work we focuse on 
the normal state temperature variation of the spectral function for the Fermions at wave vector hp as well as on the 
imaginary part of the self energy as a function of frequency. Since we are only interested here in qualitative features 
such as the temperature variation of the pseudogap, we carry out our study for the case of a ID bare density of states, 
knowing from past experience [ p2|p^ that the dimensionality only quantitavely influences the pseudogap features in 
this model. We therefore adopt the same procedure as outlined in ref |£2l but with a much denser mesh of 1000 k 
vectors in the Brillouin zone. This is required in order to track the temperature effects on the spectral function and 
to dertermine lip as a function of temperature, given the fact that for our two component system the number of 
Fermions and Bosons is not conserved separately. 

Before reporting on these numerical results, let us consider some basic features of the B-F model from which the 
underlying physics can be undestood. The self energies for the Fermions I]i?(k, iwn) and Bosons ^^(q, itJm) within 
lowest order perturbation theory are given by 

EHk,^c.„)^^^ "^^^--;^) + '-^("^°) (2) 

N '^^ lUJn- Eo+ £-k+q 

V (r, ^,, \-'"^ Sr l-»F(£-k+q)-riF(gk) , . 

i.B(q, «Wm) - T7 >, ■■ ('^) 

N ^ UjJ„i - £-k+q - £k 

where np^B{---) denote the Fermi and respectively Bose distribution function, e^ = t{z — ^^ exp{i]<.d) — fj.) is the free 
Fermion dispersion with S being the vectors linking nearest neighbour sites, Eq = Ag — 2/i and N is the number of 
k vectors in the Brillouin zone. Considering the problem of the B-F model within the framework of a selfconsistent 
RPA scheme requires to look for distribution functions n^ and n^ which replace np{...) and nB{---) in Eqs.(2,3) and 
which have to be determined selfconsistently via the equations 

with the Fermi and Bose Green's functions given by 

Gp{\<i,iujn) = [iw„ - Ck- SF(k,iw„)]"\ (5) 

Gs(q, ii^ni) = [i>^ni - Eq- SB(q, iuj,n)]~^ (6) 



Contrary to the conserving diagramatic approach to this problem which we have solved previously |2^,^, we are 
unable to find numerically stable solutions for the set of equations corresponding to the above set of selfconsistent 
RPA equations Eqs.(2-6) in 2D as well as in 3D. That is a first indication that the Fermi liquid properties in the B-F 
model are destroyed. 

Let us next attempt to understand this behaviour by closer examining those expressions for the self energies. 
Considering to begin with the situation of Eq 3> fcsT for which wc obtain 






"-F(eq) 



Eq + Eq 



(7) 



which is independent on the wavevector. Introducing the imaginary part of this contribution to the retarded self 
energy i.e., T(lu) = —2I'm'Ep{uj) we find 



T{uj) =27rw2_^„F(eq)'5(w-^o + £q) 
q 
— 2TTv'^np{Eo - u})po{Eo + ^ - uj) 



(8) 



where we have taken the limit icun -^ lu + iO. Let us now consider for simplicity a constant model density of states 
such as po{uj) = l/2zt for —zt < ui < +zt. Upon Hilbert transforming the expression for T{uj) we obtain for the 
corresponing real part of the Fermion self energy 



R(uj) = In 

^ ' 2zt 



uo - Eo 



The excitation spectrum for the Fermions is then determined by the poles of the Fermion Green's function: 



(9) 



GpiK^-k) = 



^k - Ek - R{^k) + :^T{uJk) 



(10) 



where Wk denote the complex solutions which determine the real and imaginary part of the energy spectrum. As can 
be seen from the graphical procedure for determining the poles of this Green's function, see Fig.l, we generally expect 
three poles to occur. The lowest energy real pole corresponds to the bonding state, the highest energy real pole to the 
anti-bonding state and the complex pole in the intermediary energy range to the non-bonding state. Given the fact 
that Eq > 0, since the chemical potential always must lie below the Bosonic level, the distribution function rik for the 
Fermions is enterly described by the lowest energy solution of Eq.(lO) and for which r((jjk) = 0. Wc thus obtain: 



rik = / dLL}S{uj — e]i — R{Lu))np{Lu) 
ZknF(wk) = •Z'k 



(11) 



with Z]f^= [1 — dR/du!\i^^]^^ representing the spectral weight of the lowest energy pole ojk < of the Fermion Green's 
function i.e., G_F(k, uj) — Zi^/{uj — tJk + «0). The standard way of determining the Fermi wavevector for Fermi liquids 
is to use the relation for the Fermion distribution function derived from Eq.(lO), together with the condition Wk < 
which leads to Ek^ + -^(0) — 0- The relations describing the quasi-particle properties at the Fermi surface are then 
given by: 



ek^ = - 
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(12) 



The B-F model thus predicts two well separated peaks at the Fermi wave vector whose energies are given by w ~ 
and CO o± Eo + ^^ + /^-^|^ cxp -(^o + m)/w^- 

From the above Eqs. (11,12) we notice that although in general the Fermionic distribution function has a discon- 
tinuity at ki?, this discontinuity disappears in the limit Eq ^ (which corresponds to the limit T ^ Tc{= in our 



case of a ID bare density of states) and for which e^p => oo and n-^^p ^ 0. The physical reasons for this breakdown of 
Fermi hquid properties becomes obvious by inspection of the contribution proportional to hb^Eq) in the self energy 
and which for Eq <^ ksT becomes singular. Let us for that purpose suppose that the Bosons have aquired some 
itinerancy i.e., Eq = Eq + {fiqf' /2mB where niB denotes some Boson mass [gO|. In the limit £^0^0 the term 
proportional to nB{Eq) becomes singular for q^ 0. Putting nB{Eq) — nB{Eq) + Sq^^Nno — nB{Eq)) leads to 

j:{k,uj)~v^ '^ (13) 

and ultimately to a contribution to the Ferniion Green's function, given by 

GfO^^co)^ ^ ""^^^ (14) 



This is precisely the BCS Green's function with its two poles Wk = ±7k = i-^/e^ + ij^rio showing its gap in the 
excitation spectrum given by \/v^no and a Fermion distribution function nt = |[1 — £k/7k] which has no longer a 
discontinuity at k = k^^. 

III. SPECTROSCOPIC SIGNATURES OF A TWO-COMPONENT SCENARIO 

After these phase space considerations for the spectral properties of the B-F model, let us now come to the numerical 
solution of this problem. Following the procedure discussed previously 122] we derive the Fermionic spectral function 

A(k,w) = -2/mG(k,w + iO) (15) 

at \i.p. There is no unambiguous way of determining kp for an interacting system at finite temperatures. We follow 
here the method spectroscopists use, namely to determine hp as that vector for which the distribution function 

^A{k,u;)nF{u;) (16) 

has an inversion in slope. As in our previous calculations we use here as representative parameters w = 0.1 and 
Ab = 0.4 in units of the band width 8t and a total density of charge carriers n — 2np + 2nB ~ 1. The number of 
Fermions per spin is given by np = jj X^k '^k • The values of kp determined in this way are given as a function of 
temperature in Table I, where we compare them with values derived by other criteria such as n^ = 0.5 and using 
the Luttinger sum rule J2k<k ^ ~ Nnp. 

In Fig. 2 we illustrate ACkpjOj) for a set of temperatures from which we can see the evolution of the peak structure 
of this spectral function at the Fermi level. We also observe that the lowest energy peak, corresponding to the 
bonding two-particle state, shifts to energies below the Fermi energy (equal to zero in our notation) while the non- 
bonding single particle state shifts to energies above the Fermi level, thus idicating the opening of a pseudogap as the 
temperature is reduced. These findings are in qualitative agreement with the recent experiments of ARPES on the 
underdoped high Tc materials [0,p[. 

Our theoretical predictions of the evolution of the three pole structure of the Fermionic spectral function should be 
a testing ground for 

i) the two component scenario of high Tc superconductivity 

ii) the pre-existence of bound electron pairs which would be at the origine for the deviations from Fermi liquid 
properties and responsible for the opening of the pseudogap. 

ARPES and BIS experiments should be able to see those features, inspite of the relatively poor resolution available 
for these type of experiments at present. 

The existence of the three peaks in the spectral function can be infered from the solution of the B-F model in the 
atomic limit for which they are given by e± — (eq + Eo/2) ± 7 for the lowest energy bonding and highest energy 
anti-bonding two particle states respectively and by Eq for the intermediary non-bonding single particle state. It 
should be noticed that those expressions for the two particle poles with 7 — ■\/(eo — Eq/2)^ + v^ correspond to those 
in the Green's function given by Eq.(14). 

The deviations from Landau Fermi liquid behaviour are contained in the frequency dependence of the Fermionic 
self energy. We trace the real and imaginary part of this function for different temperatures in Fig. (3). As we can 
see and as we have pointed out previously, this function does not go to zero at the Fermi level when the temperature 
decreases, as would be expected to happen for a standard Fermi liquid. We also notice from Fig. (3) that the real 



part of the Fermionic self energy is indeed very similar to the model self energy discussed above and which led to the 
form given by Eq.(9). This guarantees that the essential physics of the normal state properties of the B-F model, 
as discussed above in terms of the model density of states, is basically correct and should be rather insensitive to 
dimensionality. 

Our work totally neglects correlation effects as well as such important questions as the symmetry of the pseudogap 
and the incoherent background, seen in ARPES and associated with this pseudogap in certain direction [|l^. Such 
questions are presently being considered by us and will make the topic of some future publication. 
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Fig.l Plot of the real and imaginary part of the Fermion self energy for the constant model density of states and 
the graphic solution of the poles for the Fermion Green's function. 

Fig. 2 The temperature dependence of the Fermion spectral function ACkp^uj) at the Fermi vector showing the 
opening of the pseudogap and the evolution of a three-pole structure as the temperature is decreased (a). The 
temperature evolution of the photo emission spectrum intensity (b). T = 0.007 solid line, T = 0.01 dashed line and 
T = 0.02 dotted hne. 

Fig. 3 The temperature dependence of the imaginary part (a) and the real part (b) of the Fermion self energy at 
the Fermi vector, T = 0.007 solid line, T ^ 0.01 dashed hne and T = 0.02 dotted line. 



T 


np 


kf(a) 


kF(6) 


kp{c) 


0.020 


0.276 


145 


143 


138 


0.015 


0.230 


147 


147 


146 


0.010 


0.283 


150 


148 


142 


0.007 


0.284 


151 


149 


142 


0.006 


0.285 


152 


150 


143 



TABLE I. The temperature dependence of np and k^ in units of tt/IOOI. We compare kf as determined by d^rik/rfk^ = 
(a) with kF as determined by n^p = 5 (b) and with kF determined by the Luttinger sum rule X^j^^j^ 1 = Nnp (c). 
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